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ABSTRACT 



Context. NGC346-013 is a peculiar double-lined eclipsing binary in the Small Magellanic Cloud (SMC) discovered by the VLT- 
FLAMES survey of massive stars. 
Aims. We use spectroscopic and photometric observations to investigate the physical properties and evolutionary history of NGC346- 
013. 

Methods. Spectra obtained with VLT/FLAMES are used to construct a radial velocity curve for NGC346-013 and to characterise 
the early B-type secondary. Photometry obtained with the Faulkes Telescope South is then used to derive orbital parameters, while 
O |i spectra of the secondary are compared with synthetic spectra from TLUSTY model atmospheres. 

Results. The orbital period is found to be 4.20381(12) days, with masses of 19.1 ± 1.0 and 11.9 ± O.6A70. The primary is a rapidly 
rotating (v riX = 320 ± 30km s _1 ) late-0 dwarf while the secondary, an early-B giant, displays synchronous rotation and has filled 
its Roche lobe, implying that it was originally the more massive component with recent mass transfer 'spinning up' the primary to 
near-critical rotation. Comparison with synthetic spectra finds temperatures of 34.5kK and 24.5kK for the primary and secondary 
respectively, with the nitrogen abundance of the secondary enhanced compared to baseline values for the SMC, consistent with the 
predictions of models of interacting binaries. 

Conclusions. NGC346-013 likely evolved via non-conservative mass transfer in a system with initial masses ~22+15Mq, with 
the well-constrained orbital solution and atmospheric parameters making it an excellent candidate for tailored modelling with binary 
iy-j ■ evolution codes. This system will form a cornerstone in constraining the physics of thermal timescale mass transfer, and the associated 

^v^j ' mass transfer efficiency, in massive close binary systems. 

Cg ; Key words, stars: early-type - stars: fundamental parameters - bmanes: spectroscopic - galax.es: Magellan, Clouds 

d ■ 

1—1 ! 1. Introduction iBrown et al.l kOOl). while massive secondaries spun-up to crit- 
,— ' ical rotation during mas s transfer represent candidates for 7- 
ryH . Binary systems are of critical importance for understanding the ray burst proge nitors ( Wooslev & Hegerll2006L ICantiello et al. I 
> ; formation and evolution of massive stars. They permit both the 2QQ7) p r0 p er ties derived from massive binary systems may also 
k> ' dlr6Ct determinatlon of fundamental stellar properties such as provide key constr a i nts on the for matio n and evoluti on of wider 
X; mass and radius, and, through properties such as orbital period, steUar popu i ations dRitchie et al. Il2()lt IClark et al.ll201lh . 
U eccentricity and mass ratio, offer insight into the processes of In Ms we examine thg short . period eclipsing binary 
. * , massive star formation and dynamical interactions in the natal NG C346-O130 discovered in the VLT-FLAMES survey of mas- 
cluster. Binarity will strongly influence the evolutionary path siye starf T2] fieW centred Qn NGC346i a young> massive clus . 
both components will follow, with tidal interaction affecting ro- ter located at the centre of N66 [19^ the largest 
tational mixing (| de Mink et al. || 2009|) and transfer of mass and H n region in the Small Magellanic cloud (SMC) . The region 
angular momentum in the later evolutionary stages leading to has undergone extensive star formation in recent times, with 
a distribution of spectral types that differs significantly from NGC346 hosting the i arges t population of O-type stars in the 
single-star population models flEldndge et al. || 2008|) . Ejection SMC (Massey ^ rj^ |walborn et al. 200$ lEvansetal.l 
of the primary's hydrogen mantle during close binary evolu- 2Q06l whUe a large population of massive young ste llar objects 
tion will shape both the type of supernova and the natureof (YSQs) indicate that star formation is ongoing jsimon eit al. I 
the subsequent relativistic object ( Wellstei n & Langerl 11999; 

1 RA=00:59:30.27 , <5=-72:09:09.3 (J2000); star #782 from 

Send offprint requests to: b.ritchie@open.ac.uk iMassev et al. I d 1 989T> . who list m v = 14.46, (B — V) = —0.18, and a 

* Based on observations at the European Southern Observatory, red dening-free colour index Q — {U - B) - 0.72(73 - V) = -0.79 

Paranal, Chile in programmes 171.D-0237 and 081.D-0364 2 lEvans et al. I d2005Ll2006h . 
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l2007t iGou liermis et al. I l2010h . Studies of pre-main sequence 
stars indicate a complex history, with major star formation in 
NGC 346 starting ^6Myr ago and lasting about 3Myr, while 
an older population indicates an e arlier epoch of star forma- 
tion that took place ~10M yr ago dHennekemper et al~l [2008: 
ICignoni et al.1l20ToLl20Tlh . 

The main morphological features seen i n the FLAMES spec - 
tra of NGC346-013 were summarized by lEvans et al. I (|2006), 
who noted lines of He I, C III, N III and Si III consistent with a 
spectral type of B 1-1. 5 and velocity shifts between spectra that 
spanned ~400kms _1 . However, a broad Hell A4686 absorp- 
tion line was observed moving in the opposite sense to the other 
features, indicating the presence of a companion with an ear- 
lier spectral type, with the lower amplitude of the line indicating 
that it traces the more massive object in the system: through- 
out this paper we refer to the hotter, more massive object as the 
primary and the B-type star as the secondary, acknowledging the 
problems in applying these terms to evolved, short-period binary 
systems in which mass exchange may confuse their conventional 
meaning. 

The peculiar hot, underluminous primary and the presence 
of an apparently rotationally-broadened He II line suggest that 
NGC346-013 is seen in a brief evolutionary phase shortly af- 
ter the end of thermal-timescale Roche-lobe overflow, where 
mass transfer has reversed the initial mass ratio of the system 
and 'spun-up' the mass-gainer to rapid rotation. Such a sys- 
tem is of critical importance to understandin g the efficiency 



of mass transfer in short-period binaries ([W ellstein & L anger 
1999; We llstein et al. l200lHPetrovic et al. I2005tlde Mink et al 
2007), a process that shapes both the ongoing evolution 



of the systems and the extreme evolutionary endpoints that 
can only be reached via extensive binary interaction, such 
as Wray 977 (a ~40+1 . 4M^ Bl Ia+/neut ron star binary; 
Wells tein & Langer 1 119991: iKaner et al. I 120061). 4U 1700-37 (a 
58+2.4M 06.5 Iaf+ /compact object binary; IClark et al. 1120021) 
or the magnetar CXOU J 1647 10.2-4552 16, which formed a 
neutron star despite a minim um progenitor mass approaching 
5OM (Ritch ie et al. I I2Q10I) . We have therefore undertaken a 
programme of observations aimed at deriving accurate proper- 
ties for the two components o f NGC346-013, sup plementing the 
original dataset discussed by [Evans et al. I (2006) with an addi- 
tional 14 spectra obtained with FLAMES in 2008 and a pho- 
tometric dataset acquired during 2008-2009 with the Faulkes 
Telescope South (FTS). Details of our observations are given 
in Sect. [2l while in Se ct. [3] we use TLUSTY model atmospheres 
(e.g. iHubenv & Lanz Ill995h to estimate the temperature, grav- 
ity, rotational velocity and chemical abundances of the early-B 
secondary that dominates the FLAMES spectra. In Sect. [4] we 
derive orbital parameters from the radial velocity (RV) and pho- 
tometric datasets, and discuss the evolution of NGC346-013 and 
implications for NGC 346 in Sect. [5] 



2. Observational Details 

2.1. Spectroscopy 

A total of 32 spectra of NGC346-013 were obtained in 2006 as 
part of the V LT-FLAMES surve y of massive stars field centred 
on NGC 346 (Evans et al. 2006). Observations were taken with 
the GIRAFFE spectrograph, with the five relevant wavelength 
settings used listed in Table [jfl. An integration time of 2275s 



Table 1. FLAMES setups and lines used for RV measurement. 



Setup Coverage (A) i?(A c ) Lines 

He I A4009, He I A4026 
HeiA4121, Hei A4144 



He I A4388, He I A4471 
Si in AA4553, 4568 
Si m AA4553, 4568 
N in A4642, C m A4650 
He II A4686, He I A4713 



HR02 


3854-4049 


19600 


HR03 


4033-4201 


24800 


HR04 


4188^392 


20350 


HR05 


4340-4587 


18470 


HR06 


4538-4759 


20350 




4024 4028 4032 

wavelength (A) 



4024 4028 4032 

wavelength (A) 



Fig. 1. Two-component Gaussian fit to the He I A4026 lines near 
maximum separation at MJD54663.35 {<fr = 0.70; left panel) and 
approximately half an orbit later at MJD54665.29 (cf> = 0.16; 
right panel). 



was used, with the exception of HR06 #08 (MJD=52988. 18332) 
which used an integration time of 2500s (related to execution 
of the service mode observations). The data were reduced using 
version 1 .10 of the Giraffe B ase-Line Reduction Software (gir- 
BLDRS; Blec ha et aI71l2003h . with sky spectra subtracted from 
each target using the Starlink package DIPSO iHowarth et al. I 
2003). A further 14 spectra in HR02 mode were obtained during 
a 9-day period in July 2008 as part of a programme of follow- 
up observations of the cluster (PLHunter), each with an inte- 
gration time of 2662s. These data were pipeline reduced using 
version 3.6.8 of ESOREX0, with individual spectra shifted to the 
heliocentric frame using the IRAI0 RVCORRECT and DOPCOR 
routines. A median sky spectrum was then created from the sky 
fibres and subtracted from each object spectrum. The signal-to- 
noise ratio was in excess of ~75, varying slightly depending 
on the exact conditions when individual observations were ob- 
tained. 



3 A furthe r settin g, HR14A (covering 6391-6701A) is described by 
lEvans et al. I i2006l) but is not used here, as Ha and He I A6678 are 



strongly blended with nebular emission. Observations in HR04 mode 
were not used for RV measurement due to a lack of strong metal lines 
in this region of the spectrum, but the H<5 line was used for surface 
gravity fits, as discussed in Sec. 13. II 

4 http://www.eso.org/sci/data-processing/software/cpl/esorex.html 

5 IRAF is distributed by the National Optical Astronomy 
Observatories, which are operated by the Association of Universities 
for Research in Astronomy, Inc., under cooperative agreement with the 
National Science Foundation. 
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Table 2. Journal of observations. 







GIRAFFE 


RVsec 


RVp,i 


MJD a 


Phase 6 


Setup 


(km s _1 ) 


(km s _1 ) 


52926.08623 


0.445 


HR06 


245 ±7 




52926.13820 


0.457 


HR06 


244 ± 20 


152 ± 16 


52926.16522 


0.464 


HR06 


226 ± 10 


144 ± 15 


52926.20344 


0.473 


HR06 


211 ± 13 


111 ± 17 


52926.23039 


0.479 


HR06 


201 ± 5 


144 ± 16 


52954.13589 


0.117 


HR02 


313 ± 7 


47±6 C 


52954.16291 


0.124 


HR02 


335 ±9 




52954.18986 


0.130 


HR02 


344 ± 10 


47±7 C 


52955.08957 


0.344 


HR02 


370 ±9 




52955.11657 


0.351 


HR02 


348 ± 27 


40±9 C 


52955.14352 


0.357 


HR02 


337 ± 14 




52978.11549 


0.822 


HR05 


-59 ±5 


290 ± 7 C 


52978.14244 


0.828 


HR05 


-59 ± 12 


275 ± T 


52978.16940 


0.835 


HR05 


— 46 ± 14 


266 ±6 C 


52981.04545 


0.519 


HR03 


139 ±5 




52981.07248 


0.525 


HR03 


116 =t 8 




52981.09943 


0.532 


HR03 


96 ±5 




52981.13275 


0.539 


HR03 


91 ± 10 




52981.15970 


0.546 


HR03 


69 ± 11 




52981.18665 


0.552 


HR03 


70 ± 5 




52988.11104 


0.199 


HR06 


396 ± 17 


— 10 ± 11 


52988.15506 


0.210 


HR06 


412 ± 13 


-10 ± 13 


52988.18332 


0.217 


HR06 


406 ± 15 


— 15 ± 12 


53006.04733 


0.466 


HR05 


226 ± 5 




53006.07434 


0.473 


HR05 


194 ± 8 




53006.10131 


0.479 


HR05 


202 ± 15 




54663.31918 


0.693 


HR02 


-76 ±5 


293 ± 9 C 


54663.35111 


0.701 


HR02 


-78 ±6 


304 ± ll c 


54664.29465 


0.925 


HR02 


64 ±5 




54664.32772 


0.933 


HR02 


62 ±5 




54665.29193 


0.163 


HR02 


378 ± 5 


-2±7 C 


54666.25561 


0.392 


HR02 


326 ±6 




54666.28778 


0.400 


HR02 


294 ±5 




54666.35881 


0.416 


HR02 


286 ±5 




54666.39125 


0.424 


HR02 


266 ±5 




54669.36147 


0.131 


HR02 


342 ± 13 


28±8 C 


54669.39421 


0.138 


HR02 


356 ± 10 




54670.28291 


0.350 


HR02 


391 ±9 




54671.26508 


0.584 


HR02 


42 ±8 




54672.28464 


0.826 


HR02 


-66 ±9 


294 ± 9 C 



a MJD given at the midpoint of the integration, b For superior conjunc- 
tion at MJD 52924.21553(66) and orbital period 4.20381(12) days (see 
Sect. l4.lt . c From two-component Gaussian fits to the He I lines. 



The strong Balmer-series lines could not be used for RV 
measurement due to their broad wings, lack of sharp core, and 
presence of nebular emission, and RVs were instead measured 
using the strongest available helium and (when present) metal 
absorption lines in each spectrum: the lines used for each setup 
are listed in TableQ] The He I lines appear to show only minimal 
nebular components, but are notably asymmetric in most spec- 
tra with a strong, narrow and symmetric core arising in the sec- 
ondary and a weaker contribution from the hotter primary in the 
wings; example profiles near maximum separation are shown in 
Fig-Hi RVs for the secondary were therefore measured by fitting 
profiles to just the core of the strong He I component, with the 
derived RVs in excellent agreement with other lines of C III and 
Si III that appear free from strong contamination by the primary. 
The derived RV for each spectrum is then an error-weighted av- 
erage of individual absorption lines. Unblended lines from the 
primary are only directly visible in the HR06 data, with a com- 



parison of spectra from 2003 October 14 and 2003 December 15 
(MJDs of 52926 and 52988 respectively) clearly sh owing the an- 
tiphas e motion of the He II A4686 line (see Fig. 4 of IE vans et al. I 

2006) , with the weakness of the corresponding He II A4542 ab- 
sorption line suggests a classification around 09-9.5. However, 
the HR06 spectra sample only a small portion of the orbit (</> ~0— 
0.25), and we therefore use two-component Gaussian fits to the 
the strong He I A4026 and He I A4471 lines in the HR02 and 
HR05 spectra taken near maximum separation to further con- 
strain the amplitude of the primary. The He I profiles are clearly 
double at these epochs, and satisfactory fits to co-added spectra 
can be achievecfl A full list of derived RVs for both primary and 
secondary are listed in Table [2] 

2.2. Photometry 

Photometry was obtained with the Faulkes Telescope South 
(FTS) 2m Ritchey-Chretien Cassegrain, located in Siding 
Spring, Australia, as part of a follow-up photometric programme 
to monitor the target clusters of the FLAMES survey. FTS cur- 
rently uses a Merope camera (EM03), although prior to 2009 
January the camera used was an Apogee 'Hawkcam' (EA02). 
Both cameras were coupled with an EV2 CCD42-40DD CCD 
giving a 4.7' x 4.7' field of view and producing images of 
2048 x 2048 pixels, binned 2 x 2 to give 1024 x 1024 pixels 
at 0.278" pixel -1 . The source was observed using a Bessel V 
filter, with 100s exposures. Science images were produced using 
the Faulkes automatic pipeline, which de-biases and flat-fields 
the raw images. 

A total of 104 observations of NGC346 were obtained 
through the offline queue between 2008 January 10 and 2009 
December 23, with seeing values ranging from 074 to 3'.'1. 
Aperture differential photometry of NGC346-013 and two 
nearby cluster members, MPG 781 (=Sk 80, NGC346-001) and 
MPG 811, was performed using IRAF package APPHOT. Two 
moderately bright stars lie less than 3" to the south-east of 
NGC346-013 and, at the pixel scale of the FTS data, the three 
sources could not be separated by using aperture photometry 
alone. To overcome this a large 12-pixel aperture was used to 
include all three sources, making the assumption that the other 
two stars are non-variable and therefore any variability is solely 
due to the binary. The relative contribution of the target with re- 
spect to the other two stars was then measured by obtaining point 
spread function (PSF) photometry of one of the observations un- 
der excellent seeing, using IRAF/DAOPHOT. The two stars were 
found to contribute 34% of the light at phase </> = 0.7. 

3. Quantitative spectroscopic analysis of system 

3. 1. Atmospheric parameters 

To determine effective temperatures (T e ff) and surface gravities 
(\ogg) for the orbital analysis, physical parameters for the sec- 
ondary object were estimated using the non-LTE TL USTY model 
atmosphere grid, discussed bv lDufton et al. I (I2005I) and used ex- 
tensively in the analysis of B -type Magellanic C l oud stars from 
the F LAMES survey (e.g. iHunter et al.1 120071: iTrundle et all 

2007) . Standard techniques were used to derive the physical pa- 
rameters of the secondary. The effective temperature was esti- 
mated from the Si Ill/Si IV ionization equilibrium, and fits to 

6 However, in a number of spectra sky residuals fall within the weak 
primary component of the He I line, precluding accurate deconvolution 
despite the primary and secondary being well separated; these spectra 
are not used in the analysis. 
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1.1 




4090 4100 4110 4120 4130 
Wavelength (A) 



1.1 




0.6 I , , , , , , , , I 

4320 4330 4340 4350 4360 
Wavelength (A) 

Fig. 2. Quality of model fits to the observed He, HS and H7 lines 
(upper, middle and lower panels respectively) of both stars, with 
the B-star being 1.8 times brighter than the O-star. The absorp- 
tion features in the blue flank of the He line are due to interstellar 
Ca II A3968. 



the hydrogen lines were used to constrain the surface gravity. 
The Si III triplet at 4560A wa s used to estimate the microturbu- 
lence (see lHunter et al. l2007l for a detailed discussion). Deriving 
the atmospheric parameters is an iterative process and is further 
complicated by contamination of the spectra of the secondary by 
the O-type primary, and an estimate of the continuum contam- 
ination must therefore be included. The fitting of the hydrogen 
lines for the gravity estimate is dependent on the level of con- 
tinuum contamination from the primary object, and a secondary 




3985 3990 3995 4000 4005 

Wavelength (A) 



Fig. 3. Example fits to the metal absorption lines: (Top) Si III 
triplet (bottom) N II 3995 line. 



1.8 times brighter than the primary provides the best fit to the 
observed spectra (i.e. a continuum contamination of 35%). As 
the spectra of the primary and secondary objects are well sepa- 
rated in the HR02 and HR04 observations, the He and H7 lines 
have been used to determine the surface gravity of the secondary 
object, with the best estimates being 3.40 dex and 3.50 dex from 
the two lines respectively. The individual fits to both these lines 
as well as the H<5 line where the primary and secondary are not 
well separated are shown in Fig. [2] and provide satisfactory fits 
in all cases. We derive T c ff = 24 500 K in reasonable a greement 
with the temperature scale from iTrundle et al71 (120071) for early 
B-type giants in the SMC. A microturbulence (£) of 4 km s -1 
was estimated from the Si III triplet of lines. Conservative uncer- 
tainties on these results are AT c ff = ±1500 K, A log g = ±0.2, 
and A£ = ±3kms~ 1 . 

Simultaneously to the analysis of the secondary object, the 
spectroscopic analysis of the primary object was carried out. The 
rotational velocity of the primary is much greater than the sec- 
ondary (Sec. 13.2b which precludes the observation of the metal 
lines in its spectra. As such, the He II lines were utilised for the 
effective temperature estimate of the primary object by fitting the 
observed line with rotationally broadened model spectra over a 
range of parameters, with a best fit to the He II 4542A line found 
at an effective temperature of 34 500 K. The logarithmic gravity 
of the primary was calculated to be 3.90 dex from the fit to the 
hydrogen lines. Due to the weakness of the primary spectrum, 
errors are AT ff = ±3000 K and A logy = ±0.3 respectively. 
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Table 3. Equivalent widths for the available metallic lines of the 
secondary. 



Species 


A (A) 


EWmeas. (mA) 


EWcorr. (mA) 


[X/H\ 


N II 


3995 


62 


98 


7.54 


O II 


4069 


89 


129 


8.19 


O II 


4072 


65 


99 


8.34 


O II 


4076 


75 


117 


8.35 


Oil 


4591 


51 


80 


8.21 


Oil 


4596 


48 


75 


8.23 


Oil 


4662 


42 


66 


8.04 


Mg II 


4481 


53 


34 


6.68 


Si III 


4553 


87 


142 


6.88 


Si III 


4568 


69 


108 


6.85 


Si III 


4575 


42 


66 


6.90 


Si IV 


4116 


25 


39 


6.85 



Note: Abundances (where [X/H] = log(X/H)+12) have been calculated 
by scaling the measured equivalent widths (EW mcas .) to give corrected 
widths (EW corr .) assuming that the secondary is 1.8 times brighter than 
the primary. 

Table 4. Atmospheric parameters for the secondary of NGC 346- 
013 



NGC346-013 SMC average 



Teff (K) 

wturb (km s -1 ) 

[N/H] 

[O/H] 

[Mg/H] 

[Si/H] 



24 500 ± 1500 
3.45 ±0.2 
4±3 

7.54 ±0.18 
8.23 ±0.23 
6.68 ±0.16 
6.87 ±0.30 



6.56 
8.02 
6.72 
6.77 



Note: Mean abundances (log(X /H)±12) and the uncertainties in these 
parameters are given . For comparison the mean SMC abundances from 
IHunter et ail d2007l) are also included. 



3.2. Rotational velocities 



Similar methods to those discussed in lHunter et al. I J2007) have 
again been used to estimate the rotational velocity of both the 
primary and secondary objects, viz. convolving the model spec- 
tra with rotationally-broadened profiles until a good fit to the 
observations is found. No metal lines are apparent from the pri- 
mary, so we have used the Hell line at 4686A to estimate a 
projected rotational velocity of 320 km s _1 for this objecfl A 
model profile at the estimated parameters discussed above was 
scaled to the same equivalent width as the observations in or- 
der to derive this value. However as the profile is rotationally 
dominated, errors in the underlying model spectrum should not 
be significant; nevertheless we adopt a conservative uncertainty 
of 30 kms -1 . Assuming a continuum contamination of 35% 
from the primary we derive a projected rotational velocity of 
1 10±10 km s^ 1 for the secondary object from the Si III triplet 
of lines. 

3.3. Chemical Abundances 

Using the above parameters, chemical abundances for the sec- 
ondary object have been determined (N, O, Mg and Si). The 
equivalent widths of the observed metal lines of the secondary 





0.4 0.6 
«. (P.4.203 days) 



Fig. 4. (Top) Radial velocity curve for NGC346-013. (bottom) 
Nightfall fits to the U-band lightcurve. The solid and dashed 
lines show models with and without the hotspot discussed in the 
text; residuals are plotted for the no-hotspot model. 



7 This value is in excellent agreement with the fits to the He I lines 
used to measure the RV of the primary. 



have been measured and these measured equivalent widths 
are scaled to take into consideration the continuum contam- 
ination from the primary and then the abundances are esti- 
mated. The line-by-line equivalent widths and abundance esti- 
mates are given in Table [3] with the mean values being sum- 
marized in Table |4] Additionally we give the ba se-line chem- 
ical c omposition of the SMC as determined by iHunter et al. I 
(2007) from the analysis of 14 early B-type objects. The O, Mg 
and Si abundances are in good agreement which further sup- 
ports our derived continuum contamination parameter. Although 
the nitro gen abunda n ce is ~ 1.0 de x enhanced over that of the 
SMC, IHunter et alTI (l2007l 12009-1) determined nitrogen abun- 
dances ranging from 6.5 dex to 7.6 dex and so the enhanced ni- 
trogen abundance of the secondary in NGC346-013 is not excep- 
tional. In Fig. [3] we show model fits to example metallic lines. 

The uncertainties given in Table|4]for the abundances include 
random uncertainties arising from measurement errors as well 
as the systematic uncert ainties arising from errors in the atmo- 
spheric parameters, see lHunter et al. 1 l2007h for a full discussion 
on the derivation of these uncertainties. 
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4. Orbital Parameters 

4. 1. Radial velocity curve 

An initial estimate of the orbital period was obtained using sepa- 
rate Lomb-Scargle periodograms of the radial velocity and pho- 
tometric data, with the strongest peak in both datasets found at 
0.238 cycles/day (corresponding to an orbital period of 4.202 
days), with the 2008 FLAMES dataset therefore providing near- 
contiguous coverage of two complete orbits. Taking this ini- 
tial estimate of the orbital period as a starting point, an itera- 
tive fit using the Levenberg-Marquardt algorithm was used to 
refine the period and determine the semi-amplitude and sys- 
temic velocity of the secondary; a non-zero eccentricity was 
allowed by the model, although all best fits were consistent 
with a circular orbit. The resulting solution is generally an ex- 
cellent fit to the spectroscopic results (see the top panel of 
Fig. |4|, with a period 4.20381(12) days, superior conjunc- 
tion at MJD 52924.21553(66), semi-amplitude K 2 = 250.1 ± 
2.8km s _1 , and a systemic velocity 72 = 160.9 ± 1.9km s -1 
that is very close to the local gas velocity of 162.2 ± 1.8km s _1 
found from the strong nebular Ha line. Fixing the period and 
setting the phase as 180 degrees away from the secondary then 
yields a semi-amplitude and systemic velocity for the primary of 
Ki = 156.0 ± 4.3km s _1 and 71 = 150.8 ± 3.3km s" 1 respec- 
tively. We therefore find a mass ratio q = M2/M1 = 0.62±0.02, 

(1 ± a) 2 PK 3 
Misin 3 i = v — - = 18.0 ± 0.8M Q (1) 
27rG 

and 

M 2 sm 3 i = 11.2±0.5M Q (2) 

4.2. Photometry 

The nightfall codeS3 was used to model the light curve of 
NGC346-013. Following the analysis of Sect. [TU the temper- 
ature of the secondary was set at 24500K and a model atmo- 
sphere with logy = 3.5 was used, while the primary was set 
at 34500K and \ogg = 4, and the radius was set so that the 
B-type secondary was 1.8 x more luminous than the O-star pri- 
mary. Both temperatures were allowed to vary within the error 
range of the TLUSTY model. The eccentricity was initially set 
to zero, as implied by both the RV solution and the near-contact 
light curve, while the orbital period P and mass ratio q were 
fixed at the values derived in Sect. 14.11 and a 'third light' con- 
tribution of 34% was set as discussed in Sect. 12.21 Detailed re- 
flection was included in the model, and a linear limb-darkening 
law was assumed. Errors are internal to the model and do not in- 
clude uncertainty in the 'third light' contribution, although these 
do not strongly affect the derived inclination. 

The model converges to an inclination of 78.3 ± 1.5 de- 
grees, with a 11.2i?Q secondary that has filled its Roche lobe 
and a 4.6-R© primary which is almost completely eclipsed 
at phase 0; we conservatively estimate errors on the radii of 
±O.6i?0 and ±O.5i?0 respectively. In such a configuration the 
primary eclipse is therefore deeper than secondary eclipse, de- 
spite the smaller primary being intrinsically less luminous. The 
lightcurve fit favours a slightly lower temperature ratio than 
the TLUSTY model, with a 33 500K primary and 25 000K sec- 
ondary, but poor photometric sampling of the primary eclipse 
leads to some uncertainty as to the primary temperature, and 
with errors of ±2000K and ±1000K respectively these values 

8 http://www.hs.uni-hamburg.de/DE/Ins/PerAVichmann/Nightfall 



Table 5. Summary of orbital and physical parameters of 
NGC346-013 from the analysis in Sections |4~T1 and 1421 



Parameter 


Value 




To (MJD) 


52924.21553(66) 




P (days) 


4.20381(12) 




q{= M2/M1) 


0.624 ±0.019 




a(R Q ) 


34.4 ±0.5 




i 


78.3° ± 1.5 




e 


< 0.02 






Primary 


Secondary 


Ttlljsty (K) 


34 500 ± 3000 


24 500 ± 1500 


-f NIGHTFALL (K) 


33 500 ± 2000 


25 000 ± 1000 


log g (TLUSTY) 


3.9 ±0.3 


3.45 ±0.2 


Filling factor 




0.98 ±0.02 


R(Ro) 


4.6 ±0.5 


11.2 ±0.6 


7 (km s 1 ) 


150.8 ±3.3 


160.9 ± 1.9 


^(kms^ 1 ) 


156.0 ±4.3 


250.1 ± 2.8 


Msin 3 i (M Q ) 


18.0 ±0.8 


11.2 ±0.5 


M(M ) 


19.1 ± 1.0 


11.9 ±0.6 



are consistent with the uncertainties of the fits to the FLAMES 
spectra. However, it was not possible to reproduce the clearly- 
asymmetric secondary eclipse well with the initial model, which 
also favoured a small degree of eccentricity (e ~ 0.02 ± 0.01) 
to correctly reproduce the eclipse timing even though a non- 
circular orbit is inconsistent with the putative evolutionary state. 
The uneven secondary eclipse implies a non-uniform surface 
brightness, and this was modelled by including a 'hotspot' on 
the trailing side (i.e. such that the luminosity of the secondary 
is systematically higher prior to eclipse entry, when the hotspot 
is visible, than after eclipse when the hotspot is not in the line 
of sight). When this effect was included, the resultant fit to the 
whole lightcurve was significantly improved, with a non-zero 
eccentricity no longer required to reproduce the eclipse timing. 
The best-fit hotspot has a temperature ratio of ~ 1 .5 (T ~ 38kK) 
and a longitude —30° that imply that the hotspot is not caused 
by irradiation by the hotter primary, and instead reflects a re- 
gion where the primary wind interacts with the photosphere of 
the secondarjQ Fits to the light curve both with and without the 
inclusion of a hotspot are shown in the bottom panel of Fig. [4] 
with parameters listed in Table|5] Considerable scatter is present 
in the photometry, with some data points between phases ~0.75 
and ^0.95 appearing to favour the 'no-hotspot' model. These 
data were obtained over a two-year baseline, and we therefore 
speculate that variability may be present in this feature, while 
the true interaction is undoubtedly more complex than the sim- 
ple hotspot considered here. However, we note that omission of 
the hotspot does not significantly affect either the derived incli- 
nation or the conclusion that the secondary has filled its Roche 
lobe. 

4.3. Orbital Parameters 

Combining the radial velocity and photometric models, we find 
masses for the two components of 19.1 ± I.OMq and 11.9 ± 
0.6A/ Q respectively, with the latter value consistent with the 
~B 1 spectral type of the s econ dary, e.g. the evolution ary masses 
from lHunter et al. I d2007l) and iTrundle et al. I d2007l) . From the 

9 This is ver y similar to the effect found in the 6.6-day contact binary 
Cve OB2#5 bv lLinder et ail J2009h . 
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orbital parameters we find a separation of 34.4i? Q , with a cir- 
cular orbit (e < 0.02) as expected for a Roche-lobe filling sys- 
tem. Synchronous rotation would imply v sini of 132km s -1 , 
somewhat faster than the v sini of 110km s _1 found from the 
spectral line widths of the secondary; we return to this issue in 
Sect. 15. ll Taking the photometric radius of 11.2 ± 0.6i?o and 
secondary mass of 11.9 ± 0.6M Q gives \ogg = 3.43, in ex- 
cellent agreement with the spectroscopic value of log g = 3.45 
determined in Sect. [3] The primary radius 4.6 ± 0.5i?o implies 
log g ~ 4.4 ± 0.1, somewhat higher than inferred from the 
model fit to the hydrogen line^EL while the derived radius ap- 
pears rather low compared to galactic 09 dwarfs (cf. the list 
of detached O binaries in iGiesI I2003L or V3903 Sag., a de- 
tached, short-period 07+09 binary in which the secondary has a 
near-identical mass and temperatu re to the prima ry in NGC346- 
013 but a radius ^25% higher: IVaz et al. I [19 97!) . However, the 
derived values are broadly consistent with a 20M© model at 
SMC metallicity with a zero-age main se quence radius of 5 .2R P) 
at T eff = 36 000X and \ogg = 4.3 dBrottetal.lbOllh . and 
the photometric and model atmosphere values are in acceptable 
agreement given the limitations of the photometric dataset and 
the paucity of spectral features from the O-star. 

4.4. Alternative lightcurve model 

If the radius of the primary is not constrained to ensure the 
brightness ratio implied by the model fits to the hydrogen lines, 
we find an alternative lightcurve fit at somewhat lower incli- 
nation (71° ± 2) in which the general properties of the B-type 
secondary are broadly unchanged but the O-type primary has a 
radius 12 ± 1R Q and Roche-lobe filling factor — 0.8 ± 0.1 that 
makes it slightly larger than its companion and consequently sig- 
nificantly more luminous. Revised masses for the system then 
become 21.3 ± 1.3Af and 13.3 ± 0.8M Q respectively. In this 
configuration neither component is deeply eclipsed, while inclu- 
sion of 'hotspots' on the B-type star is again required to recreate 
the profile of the secondary eclipse. However, TLUSTY model 
fits to the Balmer series lines discussed in Sect. 13. II are strongly 
inconsistent with a more luminous primary, and it was not possi- 
ble to achieve an acceptable fit to the FLAMES spectra with this 
alternative model. Consequently we do not favour this solution, 
but additional spectroscopic and photometric observations will 
be required to conclusively resolve this possible discrepancy. 

5. Discussion 

5.1. Evolutionary state 

Few clues to the nature of the primary are found in our spec- 
tra, with the majority of features within the FLAMES coverage 
originating in the ^Bl secondary. Modelling suggests a tem- 
perature of ^34.5kK and an ^09-9.5 V spectral type for the 
primary, while the high rotational velocity of ^320km s _1 de- 
rived from the broad He II lines (see Sect. l3.2b implies that it has 
been 'spun up' by mass transferred from the secondary, as tidal 
effects in close binary systems are otherwise expected to lead to 
synchronization of rotation. The secondary appears to be rotat- 
ing somewhat slower than expected for a synchronized system, 
consistent with a predicted sp i n-down of the mas s donor during 
mass transfer dLangerl Il998t lLanger et al. I l2003h . and has al- 
most filled its Roche lobe, implying that despite its current lower 

10 This discrepancy may reflect the effects of rapid rotation and a 
lower effective equatorial surface gravity on the hydrogen line profiles. 
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Fig. 5. Location of the primary and secondary relative to illus- 
trative SMC evolutionary tracks from Maeder & Meynit] (l200lh 
appropriate for the current masses of the two components. 



mass it is in a more advanced evolutionary state. The primary 
appears somewhat underluminous for its mass (see Fig. [5]), con- 
sistent with a mass-gain er that has not rejuvenated completely 
dBraun&Langerlll995l) : significant helium enrichment indica- 
tive of a more evolved state would place it close to the Wolf- 
Rayet phase, with an implied luminosity of log(L/L Q )>5.5 in- 
consistent with our observations. We note that the system does 
not feature in the X-ray cata logue from Cha ndra observations 
of the NGC 346 region by lNaz6 et all d2003l) . but the expected 
X-ray luminosities of <10 32 erg s _1 characteristic of individual 

OB stars dClark et al. Il2~008l) are too low for detection^ unless 
X-ray emission is significantly enhanced via wind interaction. 

At first sight two stars of near-equal initial mass, i.e. M; n ; ~ 
16 + 15Af©, and an orbital period of a few days would appear 
a simple model of the precursor to the current system. If the 
initial period is very short (P < 2.5 days) such a configura- 
tion may reach co ntact and merge while both stars are still o n 
the main sequence (Wellstein et al. ObOlHPetrovic et al. 1 12005). 
but for slightly longer periods the system will evolve via quasi- 
conservative case A mass transfer as the initially more mas- 
sive star overflows its Roche lobe, potentially allowing the ex- 
change of the ~ 4M needed to reach the current 19 + 12M 
mass ratio while 'spinning up' the mass gainer to rapid rotation. 
However, in such a scenario mass transfer would not halt at this 
point, but would con tinue at very high rates (> 10 -4 AjfQyr -1 ; 
Wells tein et al. l200lh until most of the hydrogen envelope of the 
mass donor is lost: the eclipsing binary RY Scuti appears to be 
possibly the only known example of this process in action, with 
ongoing mass transfer from an 8Af Q prima ry to a 30Aip com - 
panion embedded in a dense accretion disk (Smi th et al. 11201 ll) . 
The current mass ratio of NGC346-013 would therefore place it 
in the middle of fast case A transfer, but the expected duration 
of this phase is so short (~ 10 4 years) that direct observation 
is unlikely, while the extensive circumstellar material (and pos- 
sibly B[e] star morphology) indicative of rapid mass transfer is 
not observed. 

At the end of the rapid phase of case A transfer, 'slow' trans- 
fer commences, lasting for ~ lMyr at rates ~ lO^Moyr^ 1 

11 The faintest sources in iNaze et al71 d2003h have fluxes of 
~6x 10 32 erg s _1 , assuming a distance of 59 kpc for the SMC. 
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(Wel lstein et all 1200 it lLanger et al. 1 120031) . An initial 16 + 
15M system would start the slow phase with an overlumi- 
nous lo w-mass donor and a rapidly-rotating accre tor with M ~ 
25M Q dWellstein et al. ll200lblLanger et al. II2003I) . but a higher 
mass initial config uration provides a much better match: tak- 
ing model 13 of Wellst ein & Langerl (1 19991) as a template, a 
22 + 18A/q system in an inital three-day orbit will exit the fast 
phase of mass transfer with an 11.5MQ mass donor that still 
fills its Roche lobe, very similar to the B-type star in NGC346- 
013. The model is not perfect, as at the end of the fast phase 
of mass transfer both the orbital period of 7.3 days and the 
28.5A/© mass-gainer are again inconsistent with the current 
state of NGC346-013, but the model was computed assuming 
fully conservative mass transfer with no loss of mass or an- 
gular momentum from the system. However, models including 
rotation indicate that mass transfer may be significantly now- 
conservative, with only a fraction of the transferred material re- 
tained by the accretor while the remainder is ejected from the 
system dPetrovic et al. Il2005h . Assuming an efficiency /3 ~ 50% 
and a slightly greater initial mass ratio of ~ 22 + 15M , at the 
end of the fast phase system we would again observe an 11. 5M@ 
mass donor, but the accretor would have only increased in mass 
to ~ 20A/q, in good agreement with our observations, while 
the loss of angular momentum from the system might be suf- 
ficient to restrict the increase in orbital period. The B-type star 
would be expected to have nitrogen enhanced by up to 1 dex as 
CNO-pro cessed material is exposed by the loss of th e hydrogen 
envelope dPetrovic et al. 1120051 lLanger et al. 1120081) . consistent 
with the high (but not exceptional) nitrogen abundance found in 
Sect. [33] 

Such a scenario requires confirmation from binary evolu- 
tionary models tailored to match the parameters of NGC346- 
013, but represents an encouraging match to the observations 
of the system. Tidal effects are expected to slow the rotation o f 
the mass gainer within a few 10 5 years dPetrovic et al. 1 12005). 
with the still-rapid rotation of the primary therefore suggesting 
that fast transfer has ended relatively recently. Within ~ 10 6 
years the system will return to synchronous rotation, while slow 
case A transfer will finish when core hydrogen burning ends in 
the mass donor. A second phase of rapid case AB transfer then 
commences with the onset of shell burning, spinning up the O- 
type star, still on the ma in sequence, back to near-critical rota- 
tion (lLanger et al. | |2003). During this phase the mass donor will 
lose its remaining hydrogen envelope t o leave a helium core , 
while the orbital period expands greatly dWellstein et al. II2001I) . 
leaving a rapidly rotating ~ 20 — 25M Q OB (super)giant that 
appears as an apparently-single star, with the extreme mass ra- 
tio and increasingly low luminosity of the mass donor making 
detection challenging. The initially more massive star will ulti- 
mately end its life as a type Ib/c supernova, forming a neutron 
star due to the redu ced core mass resul ting from loss of the hy- 
drogen envelope dPetrovic et al. 20051) and the ea rly termina- 
tion of shell hydrogen burning (Brown et al. 11200 H) . In contrast, 
the spun-up mass gainer may then evolve under chemically ho- 
mogeneous evolution, returning to critical rotation at the end of 
core hydrogen burning and ultimately forming a rapidly-rotating 
Wolf- Rayet star that may rep resent a future 7-ray burst progeni- 
tor (cf.|CanniIEiLiD|2Q07l). 



metric surveys (e.g. OGLE-2: IUdalski et al. Ill99~7t) . and accurate 
parameters are now availabl e for significant numbers of SMC 
binar y systems (e.g. | Harries et al. 1 120031 IHilditch et al. I 120051, 
and iNorth et al. 11201 Oh , while modell ing of the IHilditch et al l 
(2005j) dataset has been carried out bv lde Mink et al. I d2007l) . In 
the context of these samples the 19.1 ±1.0 and 11.9 ± 0.6M Q 
components of NGC346-013 are unexceptional, with s everal 
more massive systems listed by IHilditch et al. I d2005l) . who 
also identify some 28 systems as being in a post-Roche-lobe 
overflow state consiste nt with the quasi-conservative models of 
Wellstein et al. (2001). However, in these samples objects near 
~ 20Mq have luminosities 0.1-0.5 dex higher than the O-type 
primary of NGC346-013, while the 320km s _1 projected rota- 
tional velocity i s also s ubstantially higher than any object listed 
by INorth et al71 d2010l) . supporting the conclusion that we are 
observing NGC346-013 in an unusual evolutionary state. 

OGLE 5202153 (P = 4.61d, 20 + 13M ; lHarries et all 
2003) appears very similar to NGC346-013, with an 09.5 pri- 
mary displaying He II absorption that is absent in the spectrum 
of the Roche-lobe filling B0. 5 III secondary. The light curve is 
also asymmetric, although in this case the luminosities of the two 
components ar e similar, a result of a larger primary radius than 
we infer here. lde Mink et al. I d2007l) model the evolution of this 
system, and find non-conservative models (f3 = 0.5) give a better 
fit than fully-conservative mass transfer. Additional parameters 
such as chemical abundances and the projected rotational veloc- 
ity of the primary that would confirm the similarity of NGC346- 
013 and OGLE 5202153 are not known. Two f urther systems, 
OGLE 5300549 (P 0lb = 1.33d, 25 + 17M Q ; IHilditch et all 
120051) a nd OGLE 9163232 (=Hodge 53-47 P nrh = 2.21d, 26 + 
16 M^: lHarries et al. 1 12003b IMorrell et aLll2003l) are identified 
bv IHilditch et al. I d2005l) as objec ts that may have rece ntly ex- 
ited the fast phase of mass transfer. IMorrell et al. I (120031) suggest 
that the 06 V and 04-5 111(f) components of OGLE 9163232 
have undergone significant mass transfer, although the current 
2.2-day orbital period would require a particularly compact ini- 
tial configuration, and the system may ultimately merge as the 
O6 V mass gainer ends core hydrogen burning, fills its Roche 
lobe, and begins reverse mass transfer while the ma ss donor 
is still on the main sequence dWellstein et al. 1 120011: see also 
Ide Minket ani2007l) . OGLE 5300549 has a B0 primary that, 
like NGC346-013, also seems somewhat underluminous for its 
mass, whi le the system has the largest temperature ratio in the 
sample of iHilditch et"aTl d2005l) . In this case th e orbital period 
is so s hort that the rotational mixing pathway of Ide Mink et al. I 
(2009) is important, and identification as a post mass-transfer 
system may be uncertain. 

5.3. Evolutionary implications for NGC346 

The systemic velocity of NGC346-013 (~155 km s _1 ) is com- 
parable with th e radial velocities of non-binary members of 
NGC346 (e.g. lEvans et al.l 120061) . while the location 2'. 34 
(equivalent to 41 pc at the distance of the S MC) from the 
04111(f) star MPG4353 in the central core dMassev et al. I 
119891) places it well within the ~7' ionized region of NGC 346 
dRelano et aini2002lFl NGC346-013 therefore appears to be a 



5.2. NGC346-013 in the context of the SMC binary 
population 

The advent of large multi-object spectrographs has allowed the 
efficient follow-up of eclipsing binaries identified in large photo- 



12 q = 00:59:04.49, S = -72:10:24.68 (J2000) 

13 The radial d i stance from the centre of the cluster was given 
by lEvans et al. I i2006l) as 1.89'; this distance was calculated 
taking coordinates for the centre of the cluster from SIMBAD 
(http://simbad.u-strasbg.fr). 
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true member of the cluster and does not appear to have been 

eje cted. 

ISabbi et al. I (120071) report extensive subclustering in 
NGC346, with NGC346-013 apparently associated with SC- 
16, a subcluster with radius 1.6pc that likely pre-dates the for- 
mation of the main NGC 346 cluster. An age of 15 ± 2.5Myr 
is reported for SC-16 by [Sabbi et aTl d2007l). consistent with 
indep endent determinati ons of 5-15Myr ( Hennekemper et al. I 
2008) and 12.5-18Myr dCignoni et al. Il2010h . but notably dis- 
crepant with ages of 3 ± lMyr for the subclusters that trace the 
major epoch of star formation in NGC 346. However, ages of 
>10Myr are incompatible with NGC346-013 being a bona fide 
member of SC-16, as even with minimum progenitor masses 
the sy stem would start t o interact at an intermediate age (~ 
8Mvr: lLaneer et al7 2008) and comp arison with SMC evolution- 
ary tracks dMaeder & Mevnetll2001l) shows that core hydrogen 
burning would have ended a t ~ lOMyr, while at the ages re- 
ported by Sab bi et al. I (120071) andlCisnoni etal.l (1201 Oh core he- 
lium burning would also have endecQ Higher progenitor masses 
would imply earlier interaction, and it therefore seems possible 
that NGC346-013 represents a chance association with SC-16, 
and is more likely to have an age compatible with the main burst 
of star formation in NGC 346. However, a young age for SC- 
16 might nevertheless be consistent with NGC346-013 being a 
member, and any conclusive determination of the likely age and 
formation scenario for NGC346-013 and its consequent implica- 
tions for the star formation history of NGC 346 must also await 
tailored modelling with binary evolution codes. 



6. Conclusions and future work 

We find current masses of 19.1 ± l.OAf© and 11.9 ± O.6M 
for the O- and B-type components of NGC346-013 respectively, 
with comparison with tailored non-LTE model atmosphere spec- 
tra favouring a more luminous B-type star and consequently a 
compact O-type primary of 4.6 ± O.5i?0. Simple conservative 
mass transfer in an initial 16+15M Q system appears unlikely, 
and we suggest that the system began life as a ^22+15M Q 
binary, with the current B-type secondary initially the more 
massive star, and has recently exited a phase of rapid, non- 
conservative mass transfer that has spun the current primary up 
to high rotation. This scenario suggests an age of <6Myr con- 
sistent with the recent burst of star formation within NGC 346. 
Within ~l-2 x 10 6 years the system will undergo a second phase 
of fast case AB mass transfer to form a low-mass helium star + 
OB (super)giant binary with a period of months. The initially 
more massive star will ultimately form a neutron star and, if 
the system is not disrupted by a supernova 'kick', a long -period 
neutron star + OB supergiant high-mass X-ray binary, while the 
mass gainer may subsequently evolve at high rotation and may 
be a possible 7-ray burst progenitor. 

Further optical observations will allow a search for evidence 
of other spectral features originating from the primary, in or- 
der to better constrain its nature; in particular ultra-violet data 
could offer insight, as the O-type star would be expected to dom- 
inate the flux in that region. High quality photometric data would 
also resolve the uncertainty as to the radius of the primary, and 
would test the robustness of the orbital parameters presented 
here. However, despite this uncertainty the general properties of 



14 While the inclusion of rotation in the evolutionary tracks results in 
core hydrogen burning ending slightly later, the expectation that a short- 
period binary will reach synchronous rotation implies that the tracks 
with rapid rotation do not apply. 



the system are well constrained, and it therefore represents an 
excellent candidate for tailored modelling with binary evolution 
codes. 
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